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ABSTRACT 



Suppressing monopoles and vortices by introducing large chemical potentials 
for them in the Wilson action for the £77(2) lattice gauge theory, we study 
the nature of the deconfmement phase transition on N% x iV T lattices for 
iV T = 4, 5, 6 and 8 and N a = 8-16. Using finite size scaling theory, we obtain 
to = ry/v = 1.93 ± 0.03 for N T = 4, in excellent agreement with universality. 
Corresponding determinations for the N T = 5 and 6 lattices are also found to 
be in very good agreement with this estimate. The critical couplings for N T = 
4, 5, 6 and 8 lattices exhibit large shifts towards the strong coupling region 
when compared with the usual Wilson action, and suggest a lot smoother 
approach to scaling. 
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1. INTRODUCTION 



Quantum field theories need regularization schemes to control divergences. 
The regularization schemes, many different types of which have been used in 
performing calculations, do not themselves affect physics in any manner, as 
they are eliminated at the end of all calculations. Long distance physics, such 
as confinement of quarks in quantum chromodynamics or determination of 
the hadronic spectrum, is conveniently studied using the lattice regulariza- 
tion. There is a lot of freedom in defining a lattice field theory. In particular, 
a variety of different choices of the lattice action correspond to the same 
quantum field theory in the continuum. While many numerical simulations 
have been performed for the Wilson action for the gauge theories, other 
choices, some motivated by the desire to find a smoother continuum limit, 
have also been used. These actions differ merely by irrelevant terms in the 
parlance of the renormalization group: in the naive classical continuum limit 
of a — > 0, they all reduce to the same Yang-Mills action and the differing 
terms are of higher order in a. 

Investigations of the deconfinement phase transition for mixed actions, 
obtained by extending the Wilson action by addition of an adjoint coupling 
term, showed |2], |3], f|, K surprising challenges to the above notion of univer- 
sality. These actions 0, |7j are 

Sbc = E {P (l - l^FUp^j +(3 A (l- -T^p)) , (1) 

and 

5v = E(^ + Pv °p) Tr ^) . (2) 

Here Up denotes the directed product of the basic link variables which 
describe the gauge fields, U^x), around an elementary plaquette P. F and A 
denote that the traces are evaluated in the fundamental and adjoint represen- 
tations respectively and the formula TyaU = \TipU\ 2 — 1 is used. o> are Z 2 
plaquette fields associated with the plaquette P and the partition function in 
the Villain case of eq.(^J) has a sum over all possible values for each of them. 
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The first term in both the equations describes the standard SU(2) Wilson ac- 
tion whereas the second term adds an adjoint SO (3) action. For zero adjoint 
coupling, i.e, for the Wilson action, several finite temperature investigations 
have shown the presence of a second order deconfinement phase transition. 
Its critical exponents have been shown § to be in very good agreement with 
those of the three dimensional Ising model, as conjectured by Svetitsky and 
Yaffe||. The verification of the universality conjecture strengthened our an- 
alytical understanding of the deconfinement phase transition which, however, 
came under a shadow of doubt by the results for the mixed actions. Follow- 
ing the deconfinement phase transition into the extended coupling plane by 
simulating these actions at finite temperature, it was found on a range of 
temporal lattice sizes for both actions ([I]) and (|2]) that: 

a] The deconfinement transition was of second order, and in agreement 
with the conjectured universality exponents, for small values of the 
adjoint coupling. It became definitely of first order for large values (3 a 
or (3 V . 

b] The deconfinement order parameter acquired a nonzero value discon- 
tinuously at the transition point for large adjoint couplings. 

c] There was no evidence of any other separate bulk transition at those 
large adjoint couplings, as expected from the results of Refs. || 0]. In 
fact, simulations on larger symmetric lattices even suggested [TO] a lack 
of a bulk phase transition at that adjoint coupling where a first order 
deconfinement transition for a lattice of temporal size four was clearly 
seen. 



Recently it was shown [11, 12] that suppression of some lattice artifacts 
such as Z 2 monopoles and vortices do restore the universality for the action 
(H|): no first order deconfinement transition was found in the entire coupling 
plane in that case. In this paper, we address this question for the action 
([!]) in the same manner and find that unlike the Villain action, one gains 
an additional bonus. The approach to scaling seems to become smoother 
than that for the original Wilson action. The organization of the paper 
is as follows: In section 2 we define the action we investigate and briefly 
recapitulate the definitions of various observables used and their scaling laws. 



3 



We present the detailed results of our simulations in the next section and the 
last section contains a brief summary of our results and their discussion. 



2. THE MODEL AND THE OBSERVABLES 





Figure 1: The phase diagram for the action (JJ), showing the first order bulk 
phase transition lines. Taken from Ref. || but with the endpoint D as 
obtained in Ref. [10|. 



Bhanot and CreutzQ found that the lattice theory defined by the ex- 
tended action of eq.(|I]) has a rich phase structure, shown in Fig. 1. Similar 
results were obtained for the Villain action (^) by Halliday and Schwimmer|7]]. 
In either case, the (3 = axis describes an SO (3) model which has a first 
order phase transition, denoted by point A in Fig. 1. At /^(or /3y) = oo, 
the theory reduces to a Z 2 lattice gauge theory with again a first order phase 
transition at (3 cnt = |ln(l + a/2) « 0.44 [jn|. For both the mixed actions, 
these first order transitions extend into the coupling plane, as shown in Fig. 
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1 by continuous lines. These lines meet at a triple point C and continue as a 
single line of first order transitions which ends at an apparently critical point 
D. The proximity of D to the (3 a — line, which defines the Wilson action, 
has commonly been held responsible for the abrupt change from the strong 
coupling region to the scaling region for the Wilson action. It has also been 
attributed as a possible reason for the dip in the non-perturbative /5-function 
obtained by Monte Carlo Renormalization Group methods. Indeed, its rela- 
tive closeness to the (3a = line for the SU (2) theory compared to the SU (3) 
theory has been thought [14 of as a possible reason for the shallower dip in 
the former case. 

The bulk transition in the Villain form of the SO (3) gauge theory is known 
[ P~5[] to be caused by a condensation of Z 2 monopoles in the strong coupling 
phase. Defining o~ c = Y\ p& d c crp for an elementary cube c, one finds that 
a c = — 1 characterizes a monopole located in c. These monopoles are absent 
in the weak coupling region, and can be suppressed at stronger couplings 
by adding a term, AmX) c (1 ~~ °c) to the action @ and setting Xm large. 
Using this extra term, Ref. demonstrated a clear merging of the second 
order deconfinement line with a first order bulk transition line for Am = 1 
for the Villain action (0). Moreover, two separate transitions were shown to 
exist on the same lattice near the merging point, thereby shedding some light 
on the paradoxes a]-c] above. While pointing to the peculiar role the bulk 
transitions play in affecting the deconfinement transitions, these results also 
suggested that Z 2 electric current loops or vortices have to be suppressed in 
restoring the universality for the mixed action fully. Defining o\ for a link 
I = (x, ji) at a point x on the lattice in the /ith direction as a product of op of 
all those plaquettes which share the link I, i.e, oi = Yl e g t ap, one finds that 
a i = — 1 signals the link to be a part of an Z 2 electric current loop. Adding 
further a term \eJ2i (1 — a i) to the action (0) in addition to the monopole 



suppression term above, Refs. [11, 12| showed that universality is restored 



in the entire coupling plane for large A#. 

While a similar mechanism is expected to work for the Bhanot-Creutz 
mixed action as well, it is clear that both the monopole and vortex 
suppression terms added above do not depend on the gauge variables Z7 M (ar) 
directly and have to be generalized suitably for addition to Sbc- One possible 
way is to define a c and o\ by replacing o~ P in them by sign(Tr^L r p). Thus 
the mixed action with chemical potentials for these monopoles and vortices 
in that case is given by 
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Sbc,s = £ i 1 ~ l^FUp^j +P A (l- ^Tr A U P 

+AmE( 1 -0 + A £ E( 1 -^)> ( 3 ) 

where a c = n pe C sign(Tr F i7p) and crj = n pe & sign(Tr F L r P ). Comparing the 
naive classical continuum limit of eq. (§) with the standard SU(2) Yang-Mills 
action, one obtains 



1 P 2/3. 



* 4 + ^' w 

Here g u is the bare coupling constant of the continuum theory. Since the 
asymptotic scaling equation for the above mixed action with additional (ir- 
relevant) couplings Xm and Xe can be easily written down in terms of g u , it 
is clear that the introduction of a non-zero Pa, Am, or Xe does not affect the 
continuum limit: the theory for each Pa, including the usual Wilson theory 
for p A = 0.0, flows to the same critical fixed point, g° = 0, in the continuum 
limit for all Am and Xe and has the same scaling behavior near the critical 
point. 

One sees that even for Pa = eq. @ corresponds to a modified Wilson 
action with different densities of the monopoles and vortices depending on the 



values of Am and Xe respectively. By analogy with the works of Refs. JTI], [12 
for the Villain action, one expects to eliminate all bulk transition lines in Fig. 
1 by setting Am and A e large. In particular, the critical point D is expected 
to be absent in that case, leading perhaps to a lot smoother transition from 
the strong coupling region to the scaling region for those values. It has to be 
stressed though that universality has to be tested afresh for eq.@ even for 
Pa = to be sure that the above naive argument about the Am and Xe terms 
being irrelevant is correct. This is what we do in the following by determining 
a critical index of the deconfinement phase transition. We then check whether 
the passage to scaling is affected by studying the deconfinement transition 
function of the temporal lattice size. 
It is interesting to note that no similar change occurs on the fundamental 
axis (Pv = 0) in the case of the Villain action since the variables <jp are 
decoupled from the gauge variables in that case and can be integrated out 



6 



exactly for any observable depending solely on U^x). However, even in that 
case the a c and <Ji could have been defined in terms of sign(Tr F Up) and 
similar results as we obtain below would be obtained. 

We studied the deconfinement phase transition on N% x iV T lattices by 
monitoring its order parameter and the corresponding susceptibility for N T = 
4, 5, 6 and 8 and N a = 8, 10, 12, 14, 15, and 16. We used the simple 
Metropolis algorithm and tuned it to have an acceptance rate ~ 40 %. The 
expectation values of the observables were recorded every 20 iterations to 
reduce the autocorrelations. Errors were determined by correcting for the 
autocorrelations and also by the jack knife method. The observables moni- 
tored were the average plaquette, P, defined as the average of TrpUp/2 over 
all independent plaquettes, and the absolute value of the average of the de- 
confinement order parameter | 16fl , L(n), over the three dimensional lattice 



spanned by n, where L is defined by 

1 N T 

L(n) = -Tr F HU (n,T) . (5) 

T=l 

Here Uo(n, r) is the timelike link at the lattice site (n, r). In order to monitor 
the nature of deconfinement and bulk phase transitions, we also define the 
susceptibilities for both \L\ and P: 

X\l\ = Nl ((L 2 )-(|L|) 2 ) (6) 
Xp = W*N T «P 2 > - <P> 2 ) . (7) 

According to the finite size scaling theory |]17|, the peak of the \L\ (or 



plaquette) susceptibility at the location of the deconfinement (or bulk) tran- 
sition should grow on N% x N T lattices like 

^f or P oc N: , (8) 

for fixed N T . For a second order transition, to = 7/^, where 7 and v char- 
acterize the growth of the \L\ (plaquette)-susceptibility and the correlation 
length near the critical temperature (coupling)on an infinite spatial lattice. 
If the phase transition were to be of first order instead, then one expects 
[|i"8ll the exponent u = 3, corresponding to the dimensionality of the space. 
In addition, of course, the average \L\ or plaquette is expected to exhibit a 
sharp, or even discontinuous, jump and the corresponding probability dis- 
tribution should show a double peak structure in case of a first order phase 
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transition. Such an analysis of the |L| -susceptibility for the Wilson action, 
where only (3 is nonzero in eq. ([3]), yielded || an exponent uj = 1.93 ± 0.03, 
in good agreement with the corresponding value (1.965 ±0.005) for the three 
dimensional Ising model, and the universality conjecture||. Universality of 
the continuum limit of lattice gauge theories predicts a similar deconfmement 
transition belonging to the same universality class as the three dimensional 
Ising model for all values of Pa, Am and A#. 

3. RESULTS OF THE SIMULATIONS 




Figure 2: Variation of the peak height of \L\ -susceptibility with spatial size 
N a for fixed Np = 4. The curve is a fit to the eq.(J|). 

In view of the fact that even for Pa = 0, the action @) differs from the 
Wilson action in a non-trivial manner for nonzero Am and A^, we concentrate 
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here on (3a = and large Am and Xe- Our aim is to check the impact of 
the suppression of monopoles and vortices, defined by the sign(Trp[/p) as 
above, on the critical exponent u on lattices with fixed N T . We then wish to 
study the scaling behaviour of the deconfinement transition by varying iV T . 
We chose Am = 1 and Ap = 5 throughout this work. Variations with respect 
to these as well as (3a should in principle be investigated although one would 
expect the results to display universality for sufficiently large Am and \e if 
the (3a = results do so. 

3.1 N T = A 



The deconfinement phase transition on N% x 4 lattices for N a = 8, 10, 
12 and 14 lattices was studied by first making short hysteresis runs on the 
smallest lattice to look for abrupt or sharp changes in both the average pla- 
quette (P) and the order parameter (\L\). In the region of strong variation of 
(\L\), longer runs were made to check whether the |L|-susceptibility exhibits 
a peak. Histogramming technique ||19|| was used to extrapolate to nearby cou- 



plings for doing this. A fresh run was made at the x\l\ peak position and 
the process repeated until the input coupling for the run was fairly close 
to the output peak position of the susceptibility. The same procedure was 
used for the bigger lattices also but by starting from the (3 C of the smaller 
lattice. No peak was found in the vanishingly small plaquette susceptibility 
throughout, suggesting a lack of any nearby bulk transition. This should be 
contrasted with the results for the Am = As = 0, which is known to ex- 
hibit a peak. Typically 100-200 thousand measurements (2-4 million Monte 
Carlo iterations) were used to estimate the magnitude of the peak height and 
the peak location for each N a . Table 1 lists our final results for all the Aro- 
used. The errors on (3 C were estimated by varying the bin size while those for 
X max were taken to be the maximum of the errors for all bin sizes. Fitting 
the peak heights in Table 1 to eq.fl^), we obtained 

u = 1.93 ±0.03 . (9) 

Fig. |2| displays the very good quality of the fit. The critical exponent u is in 
excellent agreement with the values for both the standard Wilson action || 
and the 3-dimensional Ising model quoted in previous section. Fitting the 
peak locations (3 c ,N a by the usual finite size scaling expression, 

Pc,N,= Pc,oo + B/N^ , (10) 
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where v = 0.63 is the correlation length exponent for the 3- dimensional Ising 
model and B is a constant, we obtained /3 Cj00 = 1.326±0.006, which is shifted 
by about one from the corresponding value for the Aa^ = A^ = case which 
is 2.2986 ± 0.0006 ||. Inspired by the agreement of uo above, we assumed 
universality to be true for v as well in eq. (0). However, any reasonable 
variation of v between 0.33 and 1 changes the infinite volume extrapolation 
for (3 C by a few (~ 2-3) per cent only. Thus the shift 5(3 = (3™ son - (3 C = 




02 I ' i i i i i i ' I 

1.4 1.5 1.6 1.7 1.8 1.9 2 2.1 2.2 2.3 

P 

Figure 3: The deconfinement order parameter (\L\) as a function of j3 for 
various temporal lattices iV T . For iV T = 6, the circles (triangles) correspond 
to a 12 3 (15 3 ) spatial lattice. The inset displays average plaquette (P) on all 
iV T as a function of (3. The lines are smooth extrapolations of the data to 
guide the eye. 
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0.97 appears to be dominantly due to nonzero Am and Xe, i-e, suppression 
of monopoles and vortices. 

3.2 N T = 5, 6 and 8 

In order to minimize finite spatial volume effects, we chose to work with 
N a > 2N T always, as seen in sec. 3.1. Consequently, the full 4- volume N^N T 
increased rapidly as we increased N T . This resulted in progressive shrinking 
of the coupling interval in which the histogramming technique was reliable. 
We therefore used many longer runs in the region of strong variation of (\L\) 
to obtain the susceptibility directly and used the histogramming only for the 
finer determination of the critical coupling. Fig. |3] exhibits our results for 
(\L\) and (P) (shown in the inset) as a function of (3 for N T = 5, 6 and 
8. A deconfinement phase transition is clearly visible for all of them. The 
behaviour of the order parameter for two spatial volumes 12 3 and 15 3 for N T 
= 6 also supports the existence of a transition. This is more clearly seen 
in the corresponding x\l\ determinations, shown in Fig. [|. The plaquette, 
on the other hand, describes a smooth and unique curve for all N T and N a 
values. As the inset in Fig. |^ shows, plaquette values for all these lattices 
fall on the same curve, indicating an absence of any bulk transition. Tables 
2 and 3 list the estimated maxima of for N T = 5 and 6 for two different 
spatial volumes along with the corresponding peak locations. As seen in Fig. 
U they are rather close to the input f3 at which the long runs were made. 
Using our value for uj from eq.([|), determined for N T = 4, and the peak 
height for the smaller spatial volume, the ^ max on the bigger lattice can be 
predicted. These predictions are listed in the respective tables in the last 
column and can be seen to be in very good agreement with the direct Monte 
Carlo determinations. Alternatively, one can fit eq. (H) to the peak heights in 
Tables 2 and 3 and determine to again : 

uj= 1.99 ±0.17 for N T = 5 , and 

u= 1.80 ±0.14 for N T = 6 . (11) 

Both these determinations agree with the canonical values as well as our 
own value in eq. (|]) . Not only is the universality of the deconfinement phase 
transition thus verified on three different temporal lattice sizes, but it also 
confirms that the same physical phase transition is being simulated on them, 
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2.1 2.2 2.3 



Figure 4: The |L (-susceptibility as a function of (3 on lattices with N T = 6. 
The continuous lines are extrapolations using the histogramming technique. 
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thus approaching the continuum limit of a — > in a progressive manner by 
keeping the transition temperature T c constant in physical units. 

3.3 Scaling of T c 




0.01 



Figure 5: 1/N T as a function of (3 C . The squares are from this work while 



circles are from Ref. [^TJ. The full lines depict the 2- loop asymptotic scaling 
relation of eq. (|T2"D , normalized at N T = 8 in both cases. The dashed line 
denotes eq. (|14"D , normalized the same way. 



Table 4 lists the (3 C (estimated by extrapolating to infinite volume when- 
ever possible) for all the iV T values we used. The corresponding values for 
the usual Wilson action, i.e., Am = and A^= case, are also given in Table 
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4 along with the shifts caused by switching on these two couplings. The 
shifts decrease with increasing N T but nevertheless remains sizeable even for 
the largest lattice we used. Their decrease smoothens the approach to the 
scaling limit, as we shall see below. Fig. [5] shows aT c = iV" 1 as a function 
of the corresponding critical f3 for both our simulations with suppression of 
monopoles and vortices and the standard Wilson action (without any such 
suppression). The latter are taken from the compilation of Ref. ||21|| . The 
full lines in the figure show the 2-loop asymptotic scaling relation 

where 

and bi = — — , (13) 



u 24vr 2 ' 96vr 4 
are the first two coefficients of the perturbative /3-function for the SU(2) 
Yang-Mills theory. The curve in each case was normalized to pass through 
the N T = 8 data point. The dashed line describes a 'phenomenological' scaling 



equation which is similar to the eq.fll2|) but with the exponent increased by 
a factor of two: 

One sees deviations from asymptotic scaling for both the Wilson action 
and our action with suppression of monopoles and vortices. The deviations 
for the same range of N T seem larger for our action but then one is also 
considerably deeper in the strong coupling region of the Wilson action where 
one a priori would not have even expected any scaling behaviour. As the 
agreement of our results with the dashed line of eq.(14) in Fig. ^ shows, scal- 



ing may hold in this region of (3 for the suppressed action, since the relation 
between a and /3 in this region (or g 2 ) is similar to the asymptotic scaling re- 
lation, differing only in the exponent which will cancel in dimensionless ratios 
of physical quantities. It is clear that as (3 — > oo, the difference between the 
two actions must vanish. The shifts in Table 4 do show such a trend although 
the limiting point is not reached by iV T = 8 definitely. It seems likely though 
that the trend of evenly spaced transition points for our action will continue 
and the dashed line traced by its transition points will merge with the Wilson 
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action by N T ~ 25 or so, as suggested by its approach to the data for the 
Wilson action. If this were to be so, a much smoother approach to continuum 




Figure 6: Same as Fig. § but as a function of T/T c and on additional 
lattices, as indicated. 



limit is to be expected after the suppression of monopoles and vortices. In 
particular, one expects that dimensionless ratios of physical quantities at the 
deconfmement phase transition couplings should be constant, already from 
(3 ~ 1.33, which is the transition point for the N T = 4. 

One possible interpretation of the results in Fig. [5] is that the proximity 
of the point D in Fig. [1] for the usual Wilson action causes the nontrivial 
curvature visible in the data for the Wilson action, and consequently its 
approach to the continuum limit is not so smooth. A strong suppression 
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of monopoles and vortices, as performed here, eliminates D, resulting in a 
smoother approach to scaling. Of course, for very small lattice spacings ( or 
large (3), no significant difference between the two will be seen but for sizeable 
values of the cut-off one may expect the action with suppression to exhibit 
a better and smoother approach to the continuum limit. We intend to check 
this by measuring the glueball spectrum at the critical couplings for N T = 
4-8. In the meantime, one can try to check this hypothesis by using eq. (0) 
to convert our xift) results to x{T/T c ), function of a dimensionless 

ratio for various N T . Fig. || depicts the susceptibility x as a function of 
T/T c on lattices with N T = 5, 6 and 8. Ideally one would have expected all 
susceptibility data for the same physical volume to fall on the same curve 
for different N T . These are the three lowest curves with physical volume 
8T~ 3 but with N T = 5, 6 and 8. Unfortunately, the order parameter (\L\) is 
not ultra-violet safe; it contains divergent contributions in the continuum, 
making it N T (or a)-dependent even as a function of T/T c . Consequently, 
the corresponding x's are close but not on any universal curve. On the other 
hand, an increase in physical volume to 15.6 T -3 (the 15 3 x 6 data) and 27T~ 3 
(the 15 3 x 5 data) does seem to sharpen the susceptibility peak progressively, 
as expected. Although no quantitative analysis can be done meaningfully due 
to the cut-off dependence of the order parameter itself, the results do show 
the right trend and thus support a possible scaling in the coupling region of 
these data points. 

4. SUMMARY AND DISCUSSION 

The phase diagram of the mixed action of eq. (0) in the fundamental 
and adjoint couplings, (3 and (3a, has been a crucial input in understanding 
many properties of the SU (2) and SU (3) lattice theories and their continuum 
limits. The cross-over to the scaling region from the strong coupling region, 
as well as the dip in the non-perturbative /3-function have been attributed to 
the location of the end point D of the line of bulk first order phase transition. 
In fact, even the relative shallowness of the dip for the SU(2) case compared 
to the 577(3) case is thought to be due to the closeness of the corresponding 
end point to the (3a = Wilson axis. Adding extra irrelevant terms to 
the action one obtains the modified action of eq.([j|) in which monopoles and 
vortices can be suppressed by setting the additional couplings to large values. 
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Based on the works 111], [L2| for Villain action, one expects the phase diagram 



to change completely in that case. In particular, no phase transition lines or 
their critical end point D will be there, causing a smoother transition from 
the strong coupling region to the scaling region. 

In this paper we studied the deconfmement phase transition on the fun- 
damental axis in the Pa) coupling plane but with X M = 1 and Xe = 5, i.e., 
with strong suppression of monopoles and vortices. Our finite size scaling 
analysis yielded 1.93 ± 0.03 for the critical exponent u = j/u for lattices 
with N T = 4. This value is in excellent agreement with that|| for the Wilson 
action and the three dimensional Ising model, thus verifying the naive uni- 
versality of the modified action. However, as a result of the suppression, the 
critical coupling is shifted by about unity compared to the Wilson case. Our 
results on N T = 5 and 6 also yielded similar values for u albeit with larger 
errors, confirming that the same physical transition was being studied this 
way as a function of the lattice cut-off, a. While the aT c = A" 1 was found to 
vary slower than expected from the asymptotic scaling relation fli~2l ) for N T 
= 4-8, the data did obey a similar relation with a factor of two larger expo- 
nent. A straightforward extrapolation suggests the results from the modified 
action will merge with those of Wilson action for large N T ( of about ~ 25), 
as expected in the limit of vanishing lattice spacing a. This suggests that the 
suppression makes the approach from the strong coupling side to the scaling 
side much smoother than that for the unsuppressed Wilson action, allowing 
us to simulate the theory at smaller (3. It will be interesting to see if dimen- 
sionless ratios of physical quantities such as glueball masses or string tension 
with T c are constant in the range of critical couplings explored here. Since 
the phase diagram for SU(3), and indeed for SU(N) lattice gauge theories, 
is similar and the same mechanism is expected to work for them, it will also 
be interesting to study such suppression in those theories as well. However, 
additional possibilities for topological objects may add further complications 
and may make it necessary to suppress them as well. 
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Table 1 



The values of (3 at which long simulations were performed 
on Nl x 4 lattices, (3 C and the height of the |L (-susceptibility 
peak, xfif • 










~ .max 
X\L\ 


8 


1.37 


1.366(7) 


9.34 ± 0.07 


10 


1.344 


1.360(5) 


14.34 ± 0.11 


12 


1.331 


1.345(2) 


20.44 ± 0.29 


14 


1.34 


1.343(2) 


27.48 ± 0.64 



Table 2 



Same as Table 1 but for on iV ff 3 x 5 lattices. 





P 




A .max 
X\L\ 


^predicted 


10 


1.545 


1.570(5) 


13.17 ± 0.17 




15 


1.545 


1.558(2) 


29.57 ± 0.77 


28.78 ± 0.41 



Table 3 



Same as Table 1 but for on iV ff 3 x 6 lattices. 





P 




max 
X\L\ 


max 
Xprodicted 


12 


1.75 


1.735(5) 


16.34 ± 0.45 




15 


1.70 


1.702(2) 


24.39 ± 0.41 


25.13 ± 0.70 



Table 4 



The values of (3 C at which the deconfinement phase transition 
takes place on a lattice with temporal extension N T for eq. (§) 
for Xm = 1 and = 5 (column 2) and the usual Wilson action 
(column 3), taken from Ref.||21||. The last column lists the 
shift 5(3 = /3 e Wilson - (3 C . 







^Wilson 


5(3 


4 


1.327 ± 0.007 


2.2986 ± 0.0006 


0.972±0.007 


5 


1.56 ± 0.01 


2.3726 ± 0.0045 


0.813±0.011 


6 


1.70 ± 0.01 


2.4265 ± 0.0030 


0.727±0.010 


8 


1.933 ± 0.02 


2.5115 ± 0.0040 


0.579±0.020 



